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ABSTRACT 

I describe a tensor- vector-scalar theory that reconciles the galaxy scale success of modi- 
fied Newtonian dynamics (MOND) with the cosmological scale evidence for CDM. The 
theory provides a cosmological basis for MOND in the sense that the predicted phe- 
nomenology only arises in a cosmological background. The theory contains an evolving 
effective potential, and scalar field oscillations in this potential comprise the cold dark 
matter; the de Broglie wavelength of these soft bosons, however, is sufficiently large 
that they cannot accumulate in galaxies. The theory predicts, inevitably, a constant 
anomalous acceleration in the outer solar system which, depending upon the choice of 
parameters, can be consistent with that detected by the Pioneer spacecrafts. 
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1 INTRODUCTION 

The cosmological evidence supporting the existence of a uni- 
versal pressure-less fluid (cold dark matter or CDM) appears 
to be compelling. The amplitude of the first two acoustic 
peaks in the power spectrum of the CMB anisotropics as 
observed by WMAP implies that this fluid comprises, at 
present, about 30% of the mass density of the Universe. 
Given the even more compelling evidence for an $7 = 1 uni- 
verse, most of the remainder must consist of a negative pres- 
sure "dark energy" component that may be identified with 
a positive cosmological constant (Spergel et al. 2003, Page 
et al. 2003). It is now well-known that type la supernovae 
(Garnevitch et al. 1998, Perlmutter et al. 1999) near z ~ 0.5 
are systematically dimmer by about 0.2 magnitudes than 
would be expected in an empty, "coasting" Universe. These 
observations are often cited as evidence that the Universe 
is presently dominated by this vacuum energy. But, if we 
can trust that systematic effects are well-understood, it is 
also true that the relative brightening of SNIa beyond z ~ 1 
(Tonry et al. 2003) is equally strong evidence for matter 
domination at this relatively recent epoch, again implying 
that Q.CDM ~ 0.3. Several such independent tests have lead 
to the emergence of a "concordance model" model for the 
Universe (Ostriker and Steinhardt 1995), which is now her- 
alded as a major achievement in observational cosmology 
and in which CDM is a vital component. 

Although the large scale evidence for CDM is persua- 
sive, the hypothesis is not without problems. While there 
is no shortage of candidate particles, the nature of this 
pressure-less fluid remains elusive. Particle physics theory 
beyond the standard model provides no definitive predic- 



tion, and attempts at direct laboratory detection of CDM 
particles (neutralinos, axions, etc.) have, so far, proved un- 
successful. Moreover, the evidence for CDM on the scale 
of galaxies is far from compelling, with "complicated", but 
poorly understood, astrophysical processes being invoked to 
explain various observed phenomena such as the tightness 
of the luminosity-rotation velocity relation for disk galaxies- 
the TuUy-Fisher law- and the presence of an upper limit on 
the surface brightness of galaxies- the Freeman law (Dal- 
canton et al. 1997, van den Bosch & Dalcanton 2000). The 
problem of matching the predicted rotation curves of dark 
halos emerging from cosmic N-body simulations with those 
actually observed in low surface brightness galaxies has pro- 
voked an ongoing controversy (e.g. de Blok et al 2001). 

These same phenomena, on a galaxy scale, are well- 
described by an ad hoc modification of Newtonian dynamics, 
MOND, suggested by Milgrom (1983). Here it is proposed 
that, below a critical acceleration, ao, the effective gravi- 
tational acceleration, g, approaches y/Qnao where Qn is the 
usual Newtonian acceleration. This critical acceleration is, 
within an order of magnitude, comparable to the velocity 
of light times the present value of the Hubble parameter, 
i.e., ao ~ cHo/6. This simple modification accounts for sys- 
tematic aspects of galaxy photometry and kinematics, such 
as the Freeman law and the TuUy-Fisher law, as well as 
successfully predicting the detailed form of galaxy rotation 
curves from the observed distribution of detectable baryonic 
matter. These successes have been well-cataloged (Sanders 
and McGaugh 2002) and suggest that galaxy phenomenol- 
ogy can be understood without invoking the existence of 
dark matter (McGaugh 2004a). 

MOND is also not without problems. On a phenomeno- 
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logical level, the algorithm appears unable to account for 
the full mass discrepancy in rich clusters of galaxies; it re- 
mains necessary to invoke the presence of, as yet, unde- 
tected matter in these systems (Sanders 1999, 2003; Aguirre 
et al. 2002). Moreover, the idea has long boon criticized as 
an ad hoc, empirically based hypothesis without foundation 
in deeper theory. The absence of a viable relativistic the- 
ory has placed the issues of cosmology, structure formation, 
and gravitational lensing beyond consideration in the con- 
text of MOND; that is to say, MOND cannot be confronted 
by this entire body of observed cosmic phenomena. Recently, 
Bekenstein (2004) has written down a fully covariant theory, 
TeVeS, involving dynamic vector and scalar fields in addi- 
tion to the usual tensor field of General Relativity. This is an 
important development because the theory is free of the con- 
spicuous pathologies of earlier attempts (e.g., acausal wave 
propagation, no enhanced gravitational lensing). However, 
in the theory as it now stands, there is no obvious connec- 
tion between MOND phenomenology and cosmology as is 
suggested by the near coincidence between oo and cHq. The 
free function of the theory has two discontinuous branches- 
one for cosmology and one for mass concentrations; it would 
seem difficult to follow the cosmic growth of inhomogeneities 
in the context of a theory in which cosmology is mathemat- 
ically disconnected from inhomogeneous structure. Finally, 
the theory does not address what is perhaps the greatest 
problem for MOND: the large scale evidence for CDM cited 
above. 

The construction of a cosmologically effective theory of 
MOND, and, in particular, the reconciliation of the cosmo- 
logical evidence for CDM with the galaxy-scale success of 
MOND, is the topic of the present paper. I will show here 
that this is a possibility in the context of a generalization 
of an earlier theory also due to Bekenstein (1988a)- phase 
coupling gravity or PCG- but a generalization that requires 
a vector field as in TeVeS. 

The basic idea is this: PCG is a scalar-tensor theory 
with two scalar fields, one of which couples to matter, the 
second determining the strength of that coupling. In a cos- 
mological context, the scalar field potential is an evolving 
effective potential- which can have a well-defined minimum 
(Sanders 1989). Initially, the potential evolves faster than 
the scalar field can respond, but at some point, the scalar 
falls and oscillates about the minimum. If the bare poten- 
tial has quadratic form these oscillations constitute CDM 
with, possibly, a long Compton wavelength- so long that 
these soft bosons cannot accumulate in galaxies. At the 
present epoch, the matter-coupling field provides a fifth 
force that appears below a critical value of the scalar field 
gradient as in MOND. Cosmological CDM exists as long 
wavelength excitations of the coupling-strength field. I show 
that a preferred-frame generalization of PCG reproduces the 
galaxy scale phenomenology of MOND. I give a numerical 
example of cosmological evolution which illustrates that the 
theory can reproduce the basic properties of the concordance 
model and, in particular, is consistent with the cosmologi- 
cal evidence for CDM. I further demonstrate that solar sys- 
tem phenomenology constrains the parameters of the the- 
ory and the wavelength of the bosons, but, inevitably, this 
fifth force must be present in the outer solar system- a phe- 
nomenon possibly detected as the Pioneer effect (Anderson 
et al. 1998). 



2 A BRIEF HISTORY OF IDEAS 

The basis for the suggested scalar-tensor theories of modi- 
fied dynamics is the non-relativistic, but Lagrangian-based, 
modified Poisson equation of Bekenstein and Milgrom 
(1984): 



V • [M(|V0|/ao)V0] = 47rGp 



(1) 



where fx is the function interpolating between the Newto- 
nian regime (At(a;) = 1) and the MOND regime {fi{x) = x). 
In a scalar-tensor theory, would refer only to the scalar 
field and not to the total gravitational potential. That is to 
say, in the weak field limit, such a theory would be a two- 
field theory where, in addition to 0, there is the "normal" 
Newtonian field satisfying the usual Poisson equation. 

Phase coupling gravitation (PCG) is a covariant scalar- 
tensor generalization for such modified gravity proposed by 
Bekenstein (1988a, 1988b). Here, the scalar field is complex 
= 9^"^) and the field Lagrangian is standard: 

Ls = ^lq^4>.c.^-"+q,^q" + 2V{q)] (2) 

The non-standard aspect is that only the phase couples to 
matter jointly with the gravitational tensor g^u, i.e., 

Li = Li[exp{-2ri4,)g^^...], (3) 

This leads to the scalar field equation. 



(9 0'"). 



(4) 



Here it is obvious that the scalar amplitude squared plays 
the role of the MOND interpolating function ji. 

PCG looks like a proper field theory, with the La- 
grangian possessing a self- interaction potential V{q) as well 
as a standard kinetic term. Bekenstein demonstrated that 
MOND-like phenomenology in a static background can re- 
sult if V{q) oc — g^; i.e., galaxy phenomenology requires a 
"negative sextic" potential. In general, an attraction falling 
less rapidly than 1/r^ requires V'{q) = dV/dq < over some 
range of q (Sanders 1988). 

The theory is interesting because it can possess a viable 
cosmological limit where ao is identified with r;c(^ (Sanders 
1989), and as such, becomes a candidate for a cosmologi- 
cal effective theory of MOND. In an isotropic, homogeneous 
Universe (Robertson- Walker) , the equations for the cosmic 
evolution of the scalar phase and amplitude (Einstein con- 
formal frame) are 



= 3^0j,o 



and 



q + 3hq = -V'{q)+qtP^ 



(5) 



(6) 



where ^It, is the density parameter of the baryonic matter, 
a is the universal scale factor in terms of the present scale 
factor, time (r) is in units of Ho~^ and h is the Hubble 
parameter in units of Hq. Looking at eq. 6 we see that the 

evolution of the scalar amplitude can be described in terms 
of a time-dependent "time-like" effective potential. 



Vt{q) = V{q) + 



1 fc(r) 



2 g2 



where, from eq. 5, we have 



A;(r) = Sriflbd 



(7) 



(8) 
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The obvious attraction point of such a potential, q, is a local 
minimum, which implies that the bare potential should be 
a monotonically increasing function of q, as in V{q) = Aq". 

In a Universe with inhomogeneities, the scalar phase 
and amplitude are given, in the quasi-static limit by the 
solutions to 



V • (g'V0) 



SnriGp 



V\-qV4>-^4> = V'{q)-4>W 



(9) 



(10) 



where 4> is the cosmic time derivative of the scalar phase 
given by eq. 5 with q = q, its expectation value. Thus, the 
solution for q involves a second time-dependent "space-like" 
effective potential given by 



Vs{q) = Viq) 



1 fc(r) 



2 g4 



-q 



(11) 



That is to say, cosmology adds a timc-dcpcndcnt negative 
mass term to the bare potential. Here the attraction point q' 
is at a local maximum. If certain very general conditions on 
the bare potential are satisfied, the maximum in the space- 
like potential occurs at the same value of q as the minimum 
in the timc-liko potential, q — q; i.e., the scalar amplitude 
asymptotically approaches its cosmic value far from a mass 
concentration (Sanders 1989). A cosmological theory can be 
constructed in which the total force about a mass concen- 
tration deviates from pure Newtonian below a critical value 
of the total acceleration as in MOND, but there is a return 
to 1 /r^ attraction on larger scale with a maximum possible 
mass discrepancy oi S = 2ri^/q'^. This predicted maximum 
discrepancy is not large {5 < 10) so the theory docs not 
produce extended flat rotation curves for spiral galaxies as 
observed; this only seems possible in the context of Beken- 
stein's original V{q) oc — theory. 

Apart from this phenomenological difficulty, a more fun- 
damental problem with PCG, or any scalar-tensor theory in 
which the scalar couples to matter as a conformal factor mul- 
tiplying the Einstein metric (eq. 3), is that the scalar field 
does not interact with photons. This is to say, there is no 
enhanced gravitational lensing due to the scalar field; lens- 
ing by large astronomical systems should reveal no "dark 
matter" in sharp contrast to the observations (Bekenstein 
& Sanders 1994). This led Sanders (1997) to propose a non- 
conformal relation between the physical and Einstein met- 
rics involving an additional non-dynamical vector field A'^ 
pointing in the positive time direction in the preferred cos- 
mological frame, as in the classical stratified theories (Ni 
1972). Here the metric of the physical geometry, is re- 
lated to the gravitational metric g^u as 



-277 </> 



Qpii, — 2sinh{2iq(p)A^Av. 



(12) 



This was combined with an aquadratic Lagrangian theory 
(Bekenstein & Milgrom 1984) to produce the phenomenol- 
ogy of MOND along with enhanced gravitational lensing. 

The essential problem with this aquadratic stratified 
theory is the non-dynamic vector field; such a construct vi- 
olates the principle of General Covariancc making it im- 
possible to define a conserved energy-momentum tensor 
(Jacobson & Mattingly 2001). This has been remedied by 
Bekenstein (2004) who has endowed the vector field with its 



own dynamics and demonstrated that a theory can be con- 
structed that is causal, that produces MOND phenomenol- 
ogy about mass concentration, that possesses a sensible cos- 
mological limit, and in which the relationship between the 
deflection of photons by a mass concentration and the total 
weak-field force (scalar plus tensor) is the same as in General 
Relativity. The scalar action of Bekenstein's tensor-vector- 
scalar theory, TeVeS, is essentially that of PCG in the limit 
of weak coupling (small rf) where it can be shown that the 
kinetic term for the scalar amplitude, g, can be neglected in 
the Lagrangian; i.e., 



1 2 
2^" 



+ V{q) 



(13) 



In the weak field limit this leads to a field equation equiva- 
lent to eq. 1 with /x(x) = q^ where q is given by the solution 
to qx^ — V'{q). As a toy model, Bekenstein chooses a form 
for V{q) which yields two discontinuous branches for /i(a;); 
one for negative argument that would be relevant to cos- 
mology; and one for positive x, relevant to quasi-stationary 
systems (mass inhomogeneities). 

Here I am going to put back the kinetic term for q. 
This is vital because I will identify oscillations in q with 
CDM, so its dynamics must be followed fully. The idea that 
long-wavelength scalar field oscillations may comprise dark 
matter is also not new: Press, Ryden & Spergel (1990) pro- 
posed that such "soft bosons" would form after a late phase 
transition (after recombination) leading to the development 
of density fluctuations that would be the seeds for struc- 
ture formation. The suggested potential provides a Compton 
wavelength of 30 kpc corresponding to a supposed current 
scale of large scale structures (30 Mpc). In a unified model of 
dark energy and dark matter, Sahni & Wang (1999) also pro- 
pose scalar field oscillations as dark matter and suggest that 
a large Compton wavelength could solve several of the out- 
standing observational problems of CDM halos, such as the 
apparent absence of central density cusps and the low abun- 
dance of dwarf galaxies in the Local Group. They referred to 
such dark matter as FCDM- "Frustrated Cold Dark Mat- 
ter". This was also considered by Hu, Barkana, & Gruzi- 
nov (2000) who emphasize that the de Broglie wavelength 
(the Compton wavelength extended by c/v) is the relevant 
clustering scale for such 'fuzzy' cold dark matter. In the 
present context, the de Broglie wavelength should be suffi- 
ciently long that the bosons do not accumulate in galaxies 
at all; this could be termed SFCDM or "seriously frustrated 
cold dark matter" . 

For galaxy or solar system phenomenology, wo will see 
that it is necessary to take very weak coupling {jj « 1), 
so the scalar sector of the theory is very close to that of 
TeVeS. As in the stratified theory I will not consider here the 
separate dynamics of the vector field; although, a dynamical 
vector field is a necessary requirement for general covariance, 
and it must be included properly in any final theory. It does 
appear, however, that the only consequence of the dynamical 
vector field in the weak field limit is a rescaling of the local 
gravitational constant, G, with respect to its cosmological 
value (Bekenstein 2004, Carroll & Lim 2004, Giannios 2005). 

First of all, I describe a preferred-frame generalization 
of PCG because PCG, in a cosmological context, cannot 
produce the phenomenology of MOND. 
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3 BI-SCALAR-TENSOR- VECTOR THEORY 
(BSTV) 

3.1 The scalar field equations 

The normalized vector field, dynamical or not, takes its sim- 
plest (time-like) form in the cosmological frame. Therefore a 
theory involving such a field is inevitably a preferred frame 
theory. Moreover, the unit vector, A, may be used to form 
a second scalar field invariant 



in addition to the usual invariant 



(14) 



(Sanders 1997). Obviously J 
logical frame, and 



\ij,(l>,u (15) 
(jp in the preferred cosmo- 



K = I + J 



(16) 



becomes the square of the spatial gradient of in this frame 
{K = V0 • V(/>). This means that we can make use of the 

preferred frame to separate, and manipulate separately, the 
spatial and time derivatives of at the level of the field 
Lagrangian. 

Therefore the theory that I consider will have the gen- 
eral scalar field Lagrangian, 



q'" + h(q)K- f(q)J + 2V{q)]. 



(17) 



That is to say, separate functions of q multiply the spatial 
and temporal gradients of <f) in the preferred frame. Obvi- 
ously the fields q and cj) can no longer be identified with the 
amplitude and phase of a complex scalar, but these fields 
do play the same roles as in PCG: (j) is the matter coupling 
field and q determines the strength of that coupling. The 
interaction of with particles is described by the action 



Sp — 



-SMI'' 



dx^" dx" 
dp dp 



dp 



(18) 



where the physical metric gij,„ is given by eq. 12, and p is a 
parameter along the path. 

The scalar field dynamics comes from the action 



Ss = 



IGttG 



Ls\/—gd X 



(19) 



combined with the matter action (eq. 18 summed over par- 
ticles). Taking 5Ss = with respect to variations in 0, one 

finds 

= ^^^M- K'' + (1 + 6-*"^)^^-^] (20) 



where T^i^ is the energy-momentum tensor in the physical 
frame and 

P''" = h{q)g^^ + [h{q) - f{q)]A^A\ (21) 
For an ideal fluid this becomes 

(^'"''^.M)^. = ^(P + 3p)e-^''^ (22) 

where p and p are the density and pressure actually mea- 
sured by a co-moving observer (see Bckcnstcin 2004 for a 
derivation of the source term). We notice that, unlike usual 
scalar-tensor theory, photons and other relativistic particles 
act as a source for (j). We also notice that in the preferred 



frame P*^" takes a particularly simple form with the space- 
space {i, j) and time-time (0,0) components given by 

■ = -/(?) 



(23) 



Setting 6Ss/Sq = yields the field equation for q: 

-{g'"'<l,^);^ + lh'{q)<P,a,<l>'" 

-lih'iq) - f'{q)]A^A-^,^^,, + V'{q) = (24) 
which, in the preferred frame, simplifies to 

(p^^g.^);. - ift'(g)V0 • V0 + \nqW - V'{q) = (25) 

The complete theory includes the Einstein-Hilbert ac- 
tion for the tensor field and, in principle, an action for the 

vector field, but here the separate dynamics of the vector is 
not included. Hereafter, all operations are carried out in the 
preferred frame because this is where the equations assume 
their least complicated form. 

General empirical considerations constrain the form of 
the free functions of the theory: h{q), f(q), and V(q). First 
of all, assuming homogeneity and that e^^^"^ ~ 1 (justi- 
fied below), the cosmic time derivative of in the matter- 
dominated epoch is found by integrating eq. 22: 



(26) 



Thus, in a cosmological setting {Vq = 0, V(/> = 0), the 
evolution of q (eq. 25) is determined by a time-like effective 
potential (as in PCG) 



Vt{q) 



2 f(q) 



(27) 



with fc(r) again given by eq. 8. I wish to identify dark matter 
with q oscillations in this effective potential. The fiuid repre- 
sented by such oscillations has an equation of state p = wp, 
and, for a power law bare potential, it is straightforward to 
demonstrate that w = [n — 2)/[n + 2) (Turner 1983). Iden- 
tification of the scalar field oscillations with CDM (w = 0) 
requires n = 2; therefore, I take 

1 



V{q) 



-Aq' + B. 



(28) 



Although the theory contains an evolving component of dark 
energy (quintessence), it will be seen below that inclusion 
of the bare cosmological term B is necessary to match the 
concordance model. 

The cosmological expectation value oi q (i.e., q) is given 
by the solution to 



_ i fc(r)V(g) 

dq " 2 fiqf 



0. 



(29) 



To be an attractor, this extremum in the potential must also 

be a minimum which requires that Vt"{q) > 0. This condi- 
tion, combined with eq. 29, provides a general condition on 
/(?); i-e., 



dlnU'/f) 
dln{q) 



< 1 



(30) 



evaluated at q. 

One may question the use of such a time-dependent 
effective potential with a well defined minimum because Vt 
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and q are both functions of time. But the two timescales are 
generally very different; the expectation value q varies on a 
cosmic timescale while the period of oscillation about q is 
constant and fixed by the bare potential (oc l/\/~A). As long 
as the oscillation period is much shorter than the age of the 
Universe (not true in the early Universe) the concept of a 
time-like effective potential is meaningful. 

In a Universe with inhomogeneities, but where q can be 
neglected, we may view eq. 25 as containing a space-like po- 
tential (as in PCG) which determines the spatial dependence 
of q: 



■/(?)■ 



(31) 



The extremum at q , given by the solution to dVs/dq = 0, 
occurs at the same value of q as the extremum in Vt; i.e., q = 
q' . But the extremum in this case must bo a local maximum 
in Vs{q); i.e. Vs"{q) < 0, which implies a second condition 
on f{q); i.e.. 



dlnjf) 
dln{q) 



> 1 



(32) 



evaluated at q. If conditions 30 and 32 are met, then the cos- 
mological boundary condition on q is satisfied; i.e., far from 
a mass concentration q approaches its cosmological value. 

One choice for the free functions which satisfies these 
conditions would be h{q) = q^ and f{q) = (f . This is ap- 
pealing because it would provide a cosmological realization 
of Bekcnstoin's negative scxtic effective potential in static 
PCG and hence MOND phenomenology on an extragalactic 
scale. In fact, such a model leads to a scalar force rising as 
1/r well into a galactic mass distribution; this fifth force, in 
addition to the Newtonian gravitational force, produces not 
flat but declining rotation curves in the outer parts of spiral 
galaxies. Moreover, this theory violates precise constraints 
on attraction in the solar system. These considerations 
require a slight modification of the free functions: 



and 



/(?) = 



e2(g4+^4)- 



(33) 



(34) 



Here, in addition to the coupling strength parameter tj, a 
second parameter e (« 100?;) has been introduced. 

It may appear that the presence of three free functions, 
y(g), /i(g), and /(q), characterized by three free parameters, 
j4, ?7, and e, introduces a great deal of arbitrariness into the 
theory. In fact, given the structure of the theory (eq. 17), 
the free functions must have forms similar to eqs. 28, 33, and 
34. The potential must be quadratic to assure that w = 
for the dark matter component; in the extragalactic limit, 
g < 77, it must be the case that h[(i) — > and /(g) — > (f 
to yield MOND phenomenology in the outskirts of galaxies 
(as in PCG); and in the inner solar system, it necessary that 
/i(g) — > » rf' to assure negligible scalar field effects and 
precise inverse square attraction. This is a theory strongly 
constrained by phenomenology. 



3.2 BSTV as MOND 

The quasi-static BSTV field equation for (j) in the presence 
of mass inhomogeneities (eq. 22 with 23) is: 

SttGtjp 



V • [h(q)V^] = 



(35) 



when q < e, this becomes equivalent to that of static PCG, 

eq. 9 {h{q) « q^). 

The scalar force in the low-velocity limit is given by 



(36) 



Then we find, by identification of eq. 35 with eq. 1, that the 

MOND interpolating function /i is equivalent to q^ /2rf . In 
the limit of weak coupling, where q,aq'°' may be neglected 
in the Lagrangian as in eq. 13, and where q < q < rj (the 
outskirts of galaxies) we may neglect the first and final terms 
in eq. 25. With eqs. 33 and 34, this yields 



g 

2r?2 



V12 



(37) 



That is to say, in this limit the relation between q and V0 
becomes algebraic (as in TeVeS). Eq. 37 then implies that 
the MOND acceleration parameter should be identified as 



00 _ Vl2??|0| 



(38) 



where (f) is in units of the inverse Hubble time. In this theory, 
as in PCG, there is a clear cosmological origin of ao which 
is identified with the cosmic time derivative of the matter- 
coupling scalar field, cj). The acceleration parameter evolves 
with cosmic time, and I take ao as its current value. 

Far beyond a bounded mass distribution of total mass 
M, the scalar force (eq. 36) approaches 

/. = ^ (39) 



where 



(40) 



That is to say, the total force, scalar plus Newtonian, be- 
comes 1/r^ with a ratio of the total to Newtonian force of 
1 + 5. Because of the cosmic evolution of q it is evident that 
5 also evolves; I take go and So to be the present values. 

The full solution for q and 4> about a mass concentration 
involves solving eq. 35 simultaneously with eq. 25 which, in 
this limit, may be written as 

V'g-/i'(g)V0-V0 = K'(9) 



where the space-like effective potential is 



Vs{q) = \{Ai 



(41) 



(42) 



with h{q) and /(g) given by eqs. 33 and 34. 

The basic parameters of the theory are the scalar mass- 
squared A, the coupling constant rj, and an additional pa- 
rameter e which has no effect on galactic fields so long as 
e > lOOr? > g. The procedure followed here is to take values 
of the MOND acceleration parameter and maximum mass 
discrepancy, ao and So, that arc consistent with galaxy phe- 
nomenology in the context of MOND (i.e.. So > 25 and 
ao < 0.2) and then to determine the parameters of the the- 
ory that are required by these phenomenological constraints. 
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Figure 1. The log of the Newtonian force, /„, scalar force, /s, 
and total force, ft , in units of ao plotted against the log of the 
radial distance from a point mass in units of the MOND radius 

GM/a{). This results from numerical solution of the BSTV field 
equations with i5o = 50 and ao = 0.125. 



In fact, the values of A and rj are not specified by these 
considerations but are related through the condition that 
Vsioo) = 0. Given eq. 42 this becomes 



2. qo 2 



2 /(go) dlnjf ) 
qb^ dln{q) 



(—Y 



(43) 



where, in the approximation, I have taken /(g) ~ /{^^'n'^) 
in the limit where q < rj (the extragalactic limit) and made 
use of eq. 38 and 40. Therefore, it is necessary to assume 
a value of either A or rj, and here I take r] = 3 x 10~^ to 
be consistent with weak coupling (more physical constraints 
follow below). With So = 50, and ao = 0.125c_ffo, then A = 
lo 

10 



6.9 X 10® Ho^ (this would correspond to a mass of about 



2« ev). 

With these parameters, the solution of eqs. 35 and 41 
in the presence of a point mass yields the scalar force fs/ao 
shown in Fig. 1 as a fun ction of radius, in units of the MOND 
radius r-m = \J GM/ao. Also shown are the Newtonian force 
and the total force. We notice that when ftot/ao < 1 the 
scalar force exceeds the Newtonian force, and the total force 
falls as 1/r rather than 1/r^ 

The resulting rotation curves of spherical galaxies hav- 
ing the density distribution of a Plummer sphere (Binney & 
Tremaine 1987) and total masses of lO" M© and 3 x 10^° 
M0 are shown in Fig. 2. In both cases the core radius of 
the sphere is 2 kpc. It is evident that such a theory can re- 
produce those essential properties of MOND which so nicely 
embody the overall characteristics of galaxy rotation curves: 
no Keplerian decline beyond the visible disk; an extended 
flat rotation curve with the characteristic rotational veloc- 
ity proportional to the one-fourth power of the galax;y mass 
(the TuUy-Fisher law) ; a larger mass discrepancy in systems 
with lower surface density; and a different overall form of 



Figure 2. The rotation curves of spherical galaxies (Plummer 
spheres) having total masses of 10^^ and 3 x 10^" Mq and a core 
radius of 2 kpc in both cases. The solid lines show the rotation 
curves including the scalar force and the dashed lines show the 
Newtonian rotation curves. As in Fig. 1 this results from solution 
of the BSTV field equation for So = 50 and ao = 0.125 



rotation curves in high- and low-surface density systems. A 
difference with MOND in its original form is an eventual 
Keplerian decline in rotation velocity and, consequently, the 
presence of an upper limit to the mass discrepancy of 5-1- 1, 
but this can be large. 

Clearly this is a cosmologically effective theory of 
MOND; that is, the form of rotation curves is determined by 
the effective potential eq. 42 which depends upon the cos- 
mology via 4>. But is related to f2f, and the parameters of 
the theory as given in eq. 26; with eqs. 38 and 40 and taking 
q < rj this is 



— — S2(,eOo . 
4 



(44) 



Both from considerations of primordial rmclcosynthcsis and 
the observed CMB anisotropies it is known that f2f, « 0.04 — 
0.05. Therefore, if ao < 0.2 and (5o > 25 as is required for 
galaxy phenomenology, it must be the case that e is small. 
In the above example, (5o = 50, ao = 0.125) we find that 
e « 2 x 10"^ ff ~ 0.04. 

I now demonstrate that the theory delivering this 
MOND phenomenology can also be consistent with concor- 
dance cosmology and solar system phenomenology. 



3.3 BSTV as CDM 

The full time-like effective potential (eq. 27) is given by 



Vt(q) = \Aq^ + B+'^-k{T 



,2e'(9^ + r?^) 



(45) 
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The cosmological expectation value of g, given by the con- 
dition 29, is 



-4 



2A 



l + U + 12 



(fc(r)r 



(46) 



Given the form of k{T) is evident that q decreases as the 
Universe ages. 

The actual value of q will oscillate about this expec- 
tation value and is given by the solution to eq. 25; in the 
cosmological (homogeneous) limit this becomes 



q + ihq = -Vt'{q) 



(47) 



where Vt is determined from cq. 45. The Einstein field equa- 
tion for the gravitational metric remains the same but the 
two scalar fields contribute to the energy-momentum ten- 
sor. Thus the Ftiedmann equation for the expansion of the 
Universe becomes 

1, 



' + -[r+.f(qW + 2V{q)] (48) 



where Qr is the density parameter of radiation, and I have 
assumed that iitotai = 1- 

The response time of the scalar field to the eff'ective po- 
tential is comparable to the oscillation period at the bottom 
of the time- like potential well; i.e., for small oscillations, 

(49) 



where 



Vt"{q) = A 



1 + 



dlnif{q)/f'm 



dln{q) 



Then, from eq. 43 it follows that. 



2.2 



Soil 
ao 



(50) 



(51) 



For the parameters yielding the MOND phenomenology, de- 
scribed above, (r; = 3 x 10"^ ^ = 6.9 x IC!^ Ho^) this is 
Si 2.6 X 10"^ Ho~^ (about 3.2 x 10^ years if Ho = 75). While 
the Universe is younger than r,., the scalar field cannot re- 
spond to the evolving effective potential. The scalar field q 
will remain at its initial value until r ^ Tr, will then sock the 
effective potential minimum and oscillate about that value. 
The oscillations would thus constitute dark matter, in the 
form of soft bosons, with a Compton wavelength of about 

CTr « 110 kpC. 

The evolution of a model Universe with these properties 
has been followed by numerically integrating eqs. 26, 47, and 
48. Here I have set Vlb = 0.04 which requires e = 2 x 10~^ 
to be consistent with the supposed MOND phenomenology 
(ao ~ 0.125, So — 50). In integrating eqs. 47 and 48 it 
is assumed that the scalar field oscillations remain coher- 
ent down to the present epoch. The results are shown in 
Figs. 3-5. In Fig. 3 we see the early evolution of the scalar 
field g as a function of scale factor a. This illustrates the 
point made above: for a > 10^'\ corresponding to an age of 
roughly 3 x 10~^ Hq^^ , the response time of the scalar field 
to the effective potential becomes less than the age of the 
Universe; q feels the potential and seeks its minimum. But it 
overshoots the equilibrium point and runs against the "cen- 
trifugal" barrier k{T)^ /f{q). The oscillations persist through 
the present epoch; their contribution to the present density 
of the Universe depends upon the amplitude which depends 



ft 



o 




Figure 3. A log-log plot of the scalcir field, g, as a function of scale 
factor, a, normalized to the present scale factor, for the cosmolog- 
ical model described in the text. The heavy solid line shows the 
cosmic evolution of the expectation value, q (the bottom of the 
potential well). For a < 10~^ the response time of the scalar field 
to the effective potential is longer than the age of the Universe. 



upon the initial value of the the scalar field q. Here it was 
necessary to set this at 1.2 Mpi (Planck mass) to achieve 
^CDM ~ 0.3. The time-dependent density parameters of ra- 
diation (fir), baryons (i^b), the oscillating component of the 
scalar field {ildm) and the smooth component of the scalar 
field (fix) arc shown in Fig. 4. 

In this theory there is an evolving component of the 
vacuum energy (quintessence) equal roughly to the bare 
quadratic potential V « Aq^. With A given by eq. 43 and 
making use of eq. 40, we find that V « qo^/(5o^; i.e., the 
evolving component cannot comprise a large fraction of the 
present total energy density and be consistent with MOND. 
Therefore, I have set the bare cosmological term S = 2.1 to 
force the present composition of the Universe to be that of 
the concordance model (Oa = 0.7). 

Over the course of the evolution of the model Universe, 
from a = 10'^ to o = 2, the quantity 2ri4> changes from its 
initial assumed value of zero to —1.4 x 10~^. Therefore the 
approximation exp{—2rj())) ~ 1 in determination of (f) (eq. 26) 
appears to be justified. This also implies that the effective 
value of the constant of gravity G has not varied by more 
than a 0.5% over the history of the Universe, so there is no 
contradiction with the standard Big Bang nucleosynthesis. 

However, because of the cosmic evolution of q and (j), 
those aspects of galaxy kinematics related to MOND, the 
acceleration parameter, a, and the maximum mass discrep- 
ancy, (5, will evolve with cosmic time. This is shown in Fig. 
5 where we see these quantities as a function of scale fac- 
tor. The change in the MOND critical acceleration, a, would 
be noticeable as an evolution of asymptotic rotation veloc- 
ity and consequently of the TuUy-Fisher relation for spiral 
galaxies. For example, at a redshift of z = 1, the accelera- 
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Figure 4. The running density parameter of the various con- 
stituents of the Universe as a function of log a. Shown are Qr 
(radiation), H;, (baryons), ^dm (dark matter in the form of scalar 
field oscillations) and Ha (vacuum energy density). The solid ver- 
tical line marks the present epoch where the composition matches 
that of concordance cosmology. 



Figure 5. A log- log plot of the cosmic evolution of the MOND ac- 
celeration parameter, a = aojcHo (solid line), and the maximum 
mass discrepancy, <5 (dashed line). The MOND phenomenology 
would be clearly evident from a redshift of 2 through the present 
epoch; although the asymptotically constant value of a galaxy 
rotation velocity would increase with cosmic time. 



tion parameter is roughly half its present value which im- 
plies that the asymptotic circular velocity corresponding to 
an object of a given mass would be only 85% of its present 
value. 

The Compton wavelength of the bosonic dark matter, 
in units of the Hubble length, is given by CTr (eq. 51). Taking 
ho = 0.75 this may be written 

Ac = 8.8x 10^^ kpc (52) 

The observations of the CMB anisotropies place an up- 
per limit on Ac because the scalar dark matter must form 
well before recombination (r,. < 3 x 10^''). If the soft bosons 
are to cluster on at least the scale of the third peak in the 
CMB power spectrum, then Ac < 40 kpc. 

The condition that the bosons not collect in galaxies 
places a lower limit on Ac. The minimum clustering scale is 
the de Broglie wavelength given by Ic = \c{c/v) where v is 
the velocity dispersion of a bound object {v w 300 km/s for 
a massive galaxy). If, for example, Zc > 50 kpc it must be 
the case that Ac > .05 kpc. 

These two conditions take together restrict the range of 
r); from eq. 52 we find 

7 X 10"^^ < 77(0.125/ao)(<5o/100) < 6 x 10"^ (53) 

If rj lies within these limits, then the bosons could collect 
in clusters of galaxies {Ic < 1 Mpc) while avoiding galaxies. 
This is a possible solution to the problem of the remain- 
ing mass discrepancy in clusters in the context of MOND 
(Sanders 2003). 

Below we see that solar system phenomenology places 
an even more stringent upper limit on r] and hence, via in- 



equality 53, a lower limit on the maximum mass discrepancy 
So. 



4 SOLAR SYSTEM CONSTRAINTS 

The solution of the BSTV field equations about a galaxy- 
scale (10^^ Mq), albeit spherical, mass distribution, as 
shown in Fig. 1, implies that q « 2.2?; at the approximate 
position of the sun (r = 8 kpc), corresponding to a local 
discrepancy of 5 w 0.4. In other words, in the outer solax 
system there is a scalar field acceleration that is comparable 
to the Newtonian acceleration. 

In the inner solar system, the avoidance of locally de- 
tectable preferred frame effects limits the contribution of a 
scalar field to the total force (Sanders 1997, in preparation 
2005); specifically, the acceleration due to the scalar field 
gradient must be substantially smaller than the Newtonian 
acceleration, i.e., fs < 10^''/„. Moreover, within the orbit 
of Neptune, planetary motion restricts the total weak field 
force, Newtonian plus scalar, to be quite precisely inverse- 
square, such that any spatial variation of Kepler's constant, 
GMq, is less than several parts in about 10^ (Anderson et 
al. 1995). Therefore, the substantial scalar field acceleration 
present in the Galaxy at the radius where a test particle 
would begin to feel the acceleration of the sun (about 1 pc) 
must become significantly smaller than the Newtonian force 
within the orbit of Neptune. This provides, from the follow- 
ing argument, an upper limit on the scalar field coupling 
V- 

In the limit where q < e the quasi-static field equation 
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Figure 6. A log-log plot of the Newtonian and scalar accelera- 
tions (in units of 10~* cm/s^) about a 1 Mq point mass as a func- 
tion of distance. The dashed line is the Newtonian ax;celeration 
and the dotted, solid, and long-dashed curves are the scalar accel- 



erations corresponding respectively to »7 = 5 x 10 
and 2 x 10"^^. In all cases e = 141?? 



3 X 10" 



for q about a point mass (cq. 41) may be written 

2 i2r)GM^2 1 , -2 1^.-. 
^ 1 = <1{1J:^) - o/(9)<^ (54) 



Defining y = q/qo and x = r/rm where rm = y/GM/ao = 

\J rhTs/cio {ts is the Schwarzschild radius and th is the Hub- 
ble radius) and making use of eqs. 38 and 40, this may be 
rewritten as 



"^n r„ 2 
-Vx y 



1 <5o 



aoSo rs 



6<5o 



y- 



(55) 



2 2/3x4 

Here for the final term I have assumed that q > rj {f'{q) ~ 
2q/e^) which is generally the case within the Galaxy. 

Eq. 55 may be used to estimate the radius at which a 
massive object begins to affect the variation of q. Setting 
V^y Ri (i.e., Ay = 1) we find a critical radius given by 



l<5o 



(56) 



This means that the anomalous force about the sun would 
increase as into Tc within which the mass of the sun 

would affect the radial dependence of q, leading to a non- 
inverse square force. The requirement that this deviation 
from be unnoticeable in the inner solar system places 
a lower limit on rc (for example 1000 a.u.) and therefore a 
revised upper limit on 77; i.e.. 



^ / M ^ 0.5 / 1000 o.M.x 



(57) 



Fig. 6 illustrates the numerical solution to the BSTV 
field equations (eqs. 35 and 41) about a 1 Mq object with 
go = 27? (5o = 0.4) at infinite radius (the star is embedded 



in the Galaxy). I have taken three different values of rj (2, 
3, and 5 xlO"^^). In all cases e — 141rj. The dashed line on 
the log-log plot is the Newtonian force and the long dashed, 
dotted and solid curves show, respectively, the anomalous 
scalar force for these three different assumed values of rj. 
We see that in all cases the anomalous force increases as 
1/r^ into about 1000 A.U. and then becomes constant at 
a value between 4 and 100 xlO~* cm/s^ depending upon 
the assumed value of rj. The scalar force remains constant 
until q ^ € at which point the scalar force is given by fs ~ 
(2»7^/e^)/„; i.e., at smaller radii the force is again 1/r^ but 
a factor of 2r]^/e^ below the Newtonian force. This, in fact, 
is an additional motivation for the parameter e as given in 
eq. 33. 

It would appear that a constant force in the outer solar 
system is unavoidable in the context of this theory. This is 
interesting in view of the anomalous acceleration detected in 
the outer solar system by the Pioneer spacecrafts- the "Pio- 
neer anomaly"- which is on the order of 10~^ cm/s^ (Ander- 
son et al. 1998). If the predicted anomalous force is to be less 
than the Pioneer anomaly, then 77 < 3 x 10~^^. This upper 
limit, combined with the lower limit given by the condition 
that the soft bosons should not accumulate in galaxies, in- 
equality 53, places a limit upon the maximum discrepancy 
in the present universe; i.e., 60 > 2000 if r; < 3 x lO"'^^. 
Thus the predicted gala^xy phenomenology approaches that 
of pure MOND with a very large maximum discrepancy. 

In this limiting case, corresponding to Ac ~ 50 pc, the 
response time of the scalar field would be about 160 years; 
i.e., scalar field oscillations would appear when the Universe 
approximately 200 years old or at a = 5 x 10~®, well be- 
fore decoupling. This would seem to guarantee that the soft 
bosons could cluster down to the scale of the third or fourth 
peak in the CMB angular power spectrum, although detailed 
calculations would be required to confirm consistency with 
the observed power spectrum. This limit would also assure 
that the bosons could penetrate to the cores of rich clusters 
of galaxies. 

As in all scalar-tensor theories there is a predicted time 
variation of the baryonic mass unit or, equivalently, of the 

gravitational constant. In this case, it arises because the 
source term for the gravitational metric, g^i/, contains the 
expression exp{—2ri(j)) . Therefore, the cosmological variation 
of G would be given by 



G/G = -rjcj) = 2taoHo 



(58) 



Given that the MOND acceleration parameter Oo < 0.2 and 
that e w 200??, then with the above limit on rj this would im- 
ply that G/G ~ 10~^^ Ho, well below the current limit im- 
posed by lunar laser ranging; i.e., G/G < 0.1 Ho (Williams 
et al. 2004). 



5 CONCLUSIONS 

The near numerical coincidence between the MOND accel- 
eration parameter and cHq suggests that the proper theory 
of MOND should be an effective theory- that MOND phe- 
nomenology results only in a cosmological background. If so, 
then the goal is not to derive a relativistic theory which pre- 
dicts, on the one hand, MOND phenomenology about mass 
concentrations, and on the other hand, leads to a viable 
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cosmology in a homogeneous universe. Rather, the correct 
theory may well be one in which MOND reflects the influ- 
ence of cosmology on local particle dynamics and arises only 
in a cosmological setting. 

It goes without saying that this theory is not General 
Relativity, because in the context of General Relativity local 
particle dynamics is immune to the influence of cosmology. 
This may be demonstrated via the Birkoff theorem (Will & 
Nordtvedt 1972), but ultimately comes down to the strict 
application of the Equivalence Principle (the Strong Equiv- 
alence Principle) in GR which permits no environmental in- 
fluence on local particle dynamics (apart from tides). But, 
as was first fully appreciated by Dicke (1962), this docs not 
apply to scalar-tensor theory where the cosmology-encoded 
scalar field, determined by the universal mass distribution 
and its time evolution, pervades every corner of the Uni- 
verse. This suggests that a cosmological beisis for MOND 
may be provided by scalar-tensor theory. 

However, as was also appreciated by Dicke, if scalar- 
tensor theory is to satisfy the precise experimental con- 
straints on the universality of free fall (the Weak Equivalence 
Principle), then it must couple to matter jointly with the 
Einstein metric. The simplest such coupling, via a conformal 
transformation of the Einstein metric, leads to the serious 
problem of gravitational lensing in the context of MOND- 
no enhanced deflection of photons due to the scalar field. 
This problem may be solved by postulating a non-conformal 
relation between the Einstein and physical metrics involving 
a normalized dynamical vector field (as in TeVeS), but such 
a field may also produce unobserved effects of the preferred 
cosmological frame on local particle dynamics. 

The PCG field equation, or its weak coupling limit, 
provides a mechanism by which unwanted local dynamical 
effects of the scalar or vector fields may be suppressed in 
the limit of high field gradients (e.g., in the solar system) 
where the proper theory of gravity is very nearly GR. In 
this context the outskirts of galaxies would be the transi- 
tion region between a preferred frame tensor-vector-scalar 
cosmology and a GR-dominated local dynamics that is pro- 
tected from cosmological influence. This transition would 
be observable as an acceleration dependent deviation from 
Newtonian dynamics, or MOND. 

I have shown here that a preferred frame generaliza- 
tion of PCG, i.e. a bi-scalar-tensor vector theory (BSTV), 
can be an effective theory of MOND in the sense that it 
provides a cosmological realization of Bckcnstcin's negative 
sextic potential in the context of static PCG. Moreover, it is 
a cosmologically evolving effective potential, and the MOND 
phenomenology can only result in the context of a FRW cos- 
mology. The theory is characterized by three parameters: in 
addition to the mass-squared of the bare potential. A, and 
the coupling strength parameter, r), it was necessary to add 
a third parameter, e which permits a sensible value of Q.h in 
the presence of MOND phenomenology (oq. 44) , provides a 
return to attraction in the inner solar system (Fig. 6), 
and tames the rapid evolution of G (eq. 58). 

The evolving potential also provides a mechanism 
which, with the quadratic bare potential, inevitably pro- 
duces cosmological cold dark matter in the form of scalar 
fleld oscillations. This unavoidable appearance of cosmolog- 
ical dark matter is a primary motivation for the theory. Such 
dark matter appears to be required by a variety of consider- 



ations ranging from CMB anisotropics to the observed rela- 
tive dimming and re-brightening of SNIa. 

However, if the dark matter is not to accumulate in 
galaxies, then it must be seriously frustrated; i.e., the Comp- 
ton wavelength, extended by cjvrot should be larger than 
the scale of galaxies. I have demonstrated that this implies 
a lower limit to the strength of the scalar field coupling 
to baryons (eq. 53)- a limit which also depends upon the 
maximum ratio of the scalar-to-Newtonian forces (the max- 
imum discrepancy). There is also an upper limit on the cou- 
pling parameter required by the suppression of a non-inverse 
square anomalous force in the outer solar system- that is to 
say, any non-inverse square force must be smaller than that 
implied by the Pioneer anomaly or by planetary motion. 
However, in view of the Pioneer anomaly, it is of interest 
that the theory requires, at some level, a constant anomar 
lous acceleration in the outer solar system. 

I emphasize that the present theory is not an alternative 
to Bekenstein's TeVeS; it is an example of a TeVeS theory, 
but one which explicitly involves two scalar fields; in Beken- 
stein's theory, the second field is implicit (or not explicitly 
dynamical). However, this bi-scalar-tensor- vector theory is 
incomplete. The complete covariant theory requires consid- 
eration of the full dynamics of the vector field, and, then, fol- 
lowing Bekenstein, an examination of the properties of wave 
propagation. There remains the danger of acausal propagar 
tion or instability of the background solution. I leave this for 
later consideration because, as is evident, the procedure fol- 
lowed in constructing this theory, and that of Bekenstein, is 
different from what is usually done in relativity and cosmol- 
ogy. Here, a theory is built up from the bottom by adding 
bits and pieces as required by phenomenology. This may 
seem unfamiliar or cumbersome because the usual method- 
ology in this field is more deductive; here the approach to 
the final theory is incremental. 

It may bo that the basis of MOND lies in another direc- 
tion entirely- as a modified non-local particle action (Mil- 
grom 1994). It may also be that the observational case for 
CDM has been overstated, as argued by McGaugh (2004b). 
But the weight of present observational evidence implies 
that the final theory of MOND should reproduce, or at least 
simulate, the effects of CDM upon cosmic expansion and 
upon CMB anisotropics. 

I am grateful to Jacob Bekenstein and Moti Milgrom for 
helpful criticisms and comments which had a major positive 
impact on this work. I also thank Stacy McGaugh for helpful 
comments on the manuscript. 
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